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We report a measurement of the AjJ lifetime in the exclusive decay — > J/iphP in pp collisions at 
y/s — 1.96 TeV using an integrated luminosity of 1.0 fb~^ of data collected by the CDF II detector at 
the Fermilab Tevatron. Using fully reconstructed decays, we measure r(ASJ) = 1.593 

-0.078 (stat.) ± 

0.033 (syst.) ps. This is the single most precise measurement of t(A°) and is 3.2 a higher than the 
current world average. 

PACS numbers: 13.30.Eg 14.20.Mr 



The weak decay of quarks depends on fundamental pa- 
rameters of the standard model, including elements of the 
Cabibbo-Kobayashi-Maskawa (CKM) matrix which de- 
scribe mixing between quark families. Extraction of these 
parameters from weak decays is complicated since ob- 
served quarks are not free but are confined within color- 
singlet hadrons, as described by quantum chromodynam- 
ics (QCD). An essential tool used in this extraction is the 
heavy quark expansion (HQE) technique ,1]. In HQE, 
the total decay width of a heavy hadron is expressed as 
an expansion in inverse powers of the heavy quark mass 



rriq. Lifetime ratios of 6- flavored hadrons are predicted 
to be unity through 0{l/rnb), and 0(1/to^) corrections 
are small 2%) Detailed analysis of 0(1/to^) cor- 
rections to the lifetime ratio lead to an expected value 
of ^|gu| ~ 0.94 This theoretical prediction has been 
in poor agreement with measurements for more than a 
decade US IE IE 111 - The current world average of the 
lifetime is 1.230 ± 0.074 ps, corresponding to a ratio 

of = 0.804 ± 0.049 ^. Recent calculations includ- 

ing next-to-leading-order QCD and 0(l/m^) corrections 
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lower the prediction |lCi ll 1| , making it more consistent 
with the measurements. 

In this Letter, we present a measurement of the Ajj 
lifetime, t(A°), in the fully reconstructed decay A" 
J /iphP ^ with J /ip M'^M" Eind A° pTT~ . Charge con- 
jugate modes are implied throughout. Our data sample 
consists of 1.0 fb~^ oipp collisions at -^/s — 1.96 TeV col- 
lected by the CDF II detector at the Fermilab Tevatron 
between February 2002 and February 2006. This is the 
first measurement of t{A^^) using fully reconstructed de- 
cays that is competitive with the best previous measure- 
ments, which are based on semileptonic decays. As com- 
pared to fully reconstructed decays, measurements us- 
ing partially reconstructed semileptonic decays have ad- 
ditional uncertainties due to the missing energy of the un- 
observed neutrino and the modeling of background from 
other 6-flavored baryons. 

As a check of our method, we also measure t{B^) us- 
ing a sample of B'^ — J/ipKg decays, which has larger 
yield than the A° sample. This decay channel is topo- 
logically similar to A^ J /iphP , since both and A° 
decay with large displacement from the 6-hadron-decay 
vertex. The analysis procedure for is developed using 
_B" — > J/tpKg only and checked using other &- meson de- 
cays containing a J/tp — > fJ-'^fJ.~' in the final state. The 
Aj] lifetime was not extracted until all procedures were 
established including the estimate of the systematic un- 
certainty. 

The components of the CDF II detector relevant to this 
analysis are described briefly here: a more complete de- 
scription can be found elsewhere [l2|- Charged particles 
are reconstructed using an open-cell drift chamber called 
the central outer tracker (COT) and 7 layers of sil- 
icon microstrip detectors with radii between 2.4 cm and 
28 cm 13 • These are immersed in a 1.4 T solenoidal 
magnetic field and cover the range \ri\ < 1, where rj is 
the pseudorapiditv defined as ij = — lntan(0/2) and 9 is 
the polar angle [l^. Four layers of planar drift cham- 
bers (CMU) [13 detect muons with transverse momen- 
tum pt > 1.4 GeV/c within h j < 0.6. Additional cham- 
bers and scintillators (CMX) cover 0.6 < \ri\ < 1.0 
for muons with pT > 2.0 GeV/c. 

A sample of J/ip A'^M" candidates, collected us- 
ing a dimuon trigger, is selected to begin the reconstruc- 
tion of A° and candidates. At level 1 of a three- 
level trigger system, the eXtremely Fast Tracker (XFT) 
[isj l uses COT information to fit tracks. Those tracks 
with pt > 1.5(2.0) GeV/c are extrapolated to the CMU 
(CMX) chambers and compared with the positions of 
muon-chamber tracks. The event passes level 1 if two or 
more XFT tracks are matched to muon-chamber tracks. 
Opposite-charge and opening-angle requirements are im- 
posed at level 2. At level 3, full tracking information 
is used to reconstruct J/ip l^^fJ-^ candidates. Events 
with a candidate in the mass range 2.7 to 4.0 GeV/c^ are 
accepted at level 3 and permanently recorded for further 



analysis. 

Tracks from two triggered muons are constrained to 
originate from a common vertex to make a J/ip —^ fJ-'^ fi~ 
candidate. To ensure a high-quality vertex for the life- 
time measurement, each muon track is required to have 
at least 3 axial hits in the silicon system. The recon- 
structed fi'^ fJ-~ invariant mass is required to be in the 
range 3.014 < M^+^- < 3.174 GeV/c^. This corresponds 
to approximately ±4(t of the reconstructed width, which 
is dominated by the resolution [l9| . 

We construct Kg tt+tt" and A° — pTT~ candidates 
from pairs of oppositely-charged tracks fit to a common 
vertex. Since many Kg and A° decays occur outside some 
layers of the silicon system due to their long lifetime, their 
tracks are not required to have silicon hits. We suppress 
Kg and A" cross contamination by rejecting Kg (A") 
candidates with proton-pion (dipion) invariant mass in 
the range [1.1085, 1.1235] ([0.48175, 0.5115]) GeV/c^. 

The i?° and A° candidates are reconstructed by asso- 
ciating J/ ip candidates with Kg or A° candidates in each 
event. We choose further selection requirements for our 
&-hadron samples that optimize S/y/S + B where S and 
B are the numbers of signal and background events, re- 
spectively. In the optimization procedure, S is estimated 
using a Monte Carlo simulation, while B is estimated 
using the 6-hadron invariant mass sidebands, which are 
chosen to exclude the data used in the lifetime fits to 
avoid any potential bias. 

The selection requirements resulting from the opti- 
mization are described below. We require 0.473 < A/tttt < 
0.523 GeV/c^ and pT > 1.5 GeV/c for K^ candi- 
dates. For A" candidates, we require 1.107 < Mpj^ < 
1.125 GeV/c2 and pr > 2.6 GeV/c. We require 
L^l /ajvo > 6 for K^ and LXI/cTtvo > 4 for A°, where 

L^y is defined as the distance from the J/tp vertex to 
the y°(= Kg, A°) vertex projected onto the trans- 
verse momentum vector and (Xrvo is its estimated un- 

certainty. Both B'^ and AJJ candidates are required to 
have px > 4.0 GeV/c. Finally, the of ^ 6-hadron 
kinematic fit is required to be less than 26 for 5 degrees 
of freedom. This fit finds the best 5-hadron decay ver- 
tex and momentum subject to the constraints that the 
muon tracks originate from a common vertex, the Kg 
or A*^ tracks originate from a common vertex with com- 
bined momentum pointing back in three dimensions to 
the dimuon vertex, and the invariant mass of the two 
muons is equal to the world average J/ip mass 0. The 
invariant mass distributions of the B'^ and A" candidates 
passing these requirements are shown in Fig. ^ The 
yields are N{B° J/^I^K") = 3376 ± 88 (stat.) and 
7V(A0 ^ J/V'A") = 538 ± 38 (stat.). 

The lifetime of a and A^ is determined from 
the distribution of proper decay time t given by ct = 
Liy/{0l)T = LlyC Mfc/pJ^, where L^y is the distance 
traveled by each 6-hadron candidate along the direction 
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FIG. 1: Invariant mass distribution of (a) — > J/ipKg 
and (b) — > J/iptsP candidates. The distributions with 
ct > 200 /im, where t is the proper decay time, illustrate that 
the majority of backgrounds originate from the primary in- 
teraction vertex. The distributions are fit to the sum of a 
Gaussian signal and linear background. 

of its transverse momentum and (/37)t = p\l{cMh) is 
the transverse boost, where Mf, is the world average mass 
of the h hadron .9j. Since the J/'ip vertex ocem'S at the 
same point as the h hadron decay and is well determined, 
it is used as the &-hadron decay vertex. The h hadron 
is assumed to originate from the average beamline de- 
termined on a run-by-run basis using inclusive jet data. 
The primary vertex for a given event is the x — y position 
of this beamline at the average z coordinate of the muon 
tracks at their closest approach to the beamline. 

The lifetimes are extracted using the maximum likeli- 
hood method. The likelihood function C is multivariate, 
and is constructed from the products of single variable 
probability density functions describing the distributions 
of the invariant mass m^, cti, and their respective esti- 
mated resolutions cr™ and erf*. It is given by 

N 

i=l 



where N is the number of events in the 6-hadron mass 
window, Jb is the background fraction, and V^*, T"^ , 
and 7-"" are probability density functions for ci, cr'^*, and 
mass, respectively. The mass resolution probability dis- 
tributions T"^'" [(7^) do not appear in C because they 
are equal for signal and background, within the avail- 
able statistics. Since this is not true for the ct resolution 
distributions, must be included in C 

The mass distribution is modeled as the sum of a 
Gaussian signal and linear background, where the Gaus- 
sian width CT™ is scaled by a parameter to account for 
misestimation of the mass resolutions. The ct distri- 
bution is modeled by the sum of five components, all 
convoluted with a Gaussian resolution function with a 
scale factor parameter for the af: a positive exponen- 
tial (e~'^*'/'^^/cr) for the signal, a (5— function represent- 
ing the zero-lifetime component, one negative and two 
positive exponentials accounting for mismeasured decay 
vertices and background from other heavy- flavor decays. 
The relative contribution of each of the background com- 
ponents is determined by the fit. The ct"^' distribution is 
modeled by a Gaussian convoluted with an exponential 
for both signal and background. 

We fit over the mass range [5.170, 5.390] and 
[5.521, 5.721] GeV/c2 for and A^, respectively. These 
ranges provide a sufficient sideband to constrain the 
background shape while avoiding regions where the mass 
distribution has complex structure. For both and 
A", we require ct™ < 20 MeV/c^ and fit over the range 
[-2000, 4000] Aim in ct and [0, 100] ^im in ct=*. We 
maximize the likelihood to determine the best values of 
all fit model parameters, including the signal lifetimes 
ct(B°) = 456.8 and cr(AO) = 477.6 ^23:° Mm. 

Fit projections are shown in Fig. [21 

Systematic uncertainties come from four main sources: 
fitting procedure and model, primary vertex determina- 
tion, alignment of detector elements, and Kg or A^ point- 
ing requirement in the B'^ or A" kinematic fit. The fitting 
bias is determined using a simple Monte Carlo simula- 
tion in which events are distributed according to the fit 
model. The bias was found to be less than 0.4 /im and 
0.5 ^m for and A°, respectively. The systematic un- 
certainties due to ct resolution and mass resolution are 
estimated by including additional Gaussian components 
to their respective parts of the model in separate fits to 
the data and observing the deviations from the nominal 
result. We estimate the contribution from our mass back- 
ground model by fitting with a uniform rather than lin- 
ear background shape. The systematic uncertainty due 
to the ct background model is estimated by fitting with 
two or four background exponentials convoluted with the 
resolution function and fitting with two, three, or four 
background exponentials without convolution. We study 
a possible mass dependence in the ct background shape 
by separately fitting for B^ (A") lifetime in the following 
low and high mass regions: [5.170, 5.3225] ([5.521, 5.651]) 
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and [5.2375, 5.390] ([5.591, 5.721]) GeV/c^. The ob- 
served shifts are consistent with the statistical differences 
of the two samples for both modes. We use the average 
shift of 1.9 /xm and 4.1 /im for and A", respectively, as 
an estimate of the systematic uncertainty due to a mass- 
dependent ct background. We estimate the systematic 
uncertainty due to our u'^* and ct™ distribution models 
by the observed shift between simple Monte Carlo events 
generated with the data distributions but fit with our 
model compared with simple Monte Carlo events both 
generated and fit with our model. Using the same simple 
Monte Carlo technique, we estimate the systematic un- 
certainty due to a correlation between the ct and ct"^* for 
the background component by generating simple Monte 
Carlo data sets with the correlation observed in the data 
and fitting with our baseline model where this correlation 
is absent. 

We estimate the systematic uncertainty due to our 
primary vertex determination by comparing different 
choices of the z coordinate used to evaluate the run- 
averaged beamline. We estimate uncertainties due to 
any residual misalignments of the silicon detector us- 
ing Monte Carlo samples generated with radial displace- 
ments of individual sensors (internal alignment) and rel- 
ative translation and rotation of the silicon detector with 
respect to the COT (global alignment). We also study 
the resolution and bias on the V'^ pointing to the J/tp 
vertex in data. If these were strongly ct-dependent, the 
kinematic fit quality requirement could bias the 6-hadron 
lifetime. We observe no lifetime bias and assign uncer- 
tainties of 0.6 /im for 5° and 5.4 /im for A° based on the 
statistical precision of our study. 

The systematic uncertainties are summarized in Ta- 
ble m We obtain total systematic uncertainties of 4.9 
/im for i?" and 9.9 /im for by adding the individual 
uncertainties in quadrature. 

A number of cross checks on the analysis procedure 
are performed. We measure the and lifetimes 
which are statistically consistent with the world aver- 



TABLE L Systematic uncertainties (in /im) for the measure- 
ment of ct{B^) and cr(Ab). The total uncertainties are the 
individual uncertainties added in quadrature. 



Source 


ct{B°) 


cr(A«) 


Fitter Bias 


0.4 


0.5 


Fit Model: 






ct Resolution 


3.1 


5.5 


Mass Signal 


0.7 


2.3 


Mass Background 


0.1 


0.1 


ct Background 


0.5 


0.7 


a"* Distribution Modeling 


0.1 


0.2 


ct'" Distribution Modeling 


0.6 


0.2 


Mass-ct Background Correlation 


1.9 


4.1 


ct-a'^^ Background Correlation 


0.3 


1.3 


Primary Vertex Determination 


0.2 


0.3 


Alignment 






SVX Internal 


2.0 


2.0 


SVX/COT Global 


2.2 


3.2 


V° Pointing 


0.6 


5.4 


Total 


4.9 


9.9 



age values in the following decay modes: 
B° ^p('^K*°{K*° K+TT-), B+ iIj^'^K+, and 
B+ i;K*+{K*+ i^g7r+), with V^'' /i+/^" and 
ip' —> ^pTT'^TT^ . We search for unexpected lifetime depen- 
dence on the y° and &-hadron kinematics, data-taking 
period, number of tracks in the event, and use of silicon 
hits on daughter and muon tracks; no dependence is 
observed. Finally, we determine the lifetime using two 
alternative techniques which give results consistent with 
our baseline fit: a ci-only binned fit applied to sideband- 
subtracted data and a fit to the mass distribution in ct 
bins which is insensitive to details of the ct background 
shape. 

In summary, we measure 

r(Ag) = 1.593 toml (stat.) ± 0.033 (syst.) ps. 

As a cross check, we also measure t{B^) = 1.524 ± 
0.030 (stat.) ± 0.016 (syst.) ps which is consistent with 
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the world average t{B^) = 1.530 ± 0.009 ps. Our mea- 
surement of t(A?) is consistent with the D0 result in the 
same channel |6j at the 1.7a level and is the first mea- 
surement using a fully reconstructed mode that reaches 
a precision comparable with the previous best measure- 
ments based on semileptonic decays of the A^. It is also 
comparable in precision to the current world average, but 
is 3.2 a higher jQi]. Forming a ratio with the world average 
5° lifetime, we determine 

= 1.041 ± 0.057 (stat. + syst.). 

This ratio is consistent with the higher end of the theory 
predictions |l[lS[l3- 
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